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Preface 


THE main purpose of this book is to give a thoroughly elemen- 
tary account of the ‘operational method’ for solving linear 
differential equations with constant coefficients, subject to 
prescribed initial conditions. 

No previous knowledge of differential equations will be 
demanded from the reader, and accordingly the basic proper- 
ties of linear differential equations with constant coefficients 
are treated in the first chapter. This chapter may be found 
useful even by those readers who do not wish to learn the more 
specialized technique of operators, treated in the second 
chapter. 

The book is addressed primarily to students of the exact 
sciences and of engineering, but it deals only with mathema- 
tical techniques and does not contain any illustrative examples 
drawn from other fields (such as mechanics, electric circuit 
theory, probability theory, and so on). Each reader will no 
doubt meet problems in his own field of studies to which the 
mathematical techniques described in this book may usefully 
be applied. 

Iam very grateful to Dr. Walter Ledermann for reading the 
draft with great care and for making many useful suggestions. 
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CHAPTER ONE 


Linear Differential Equations with Constant 
Coefficients 


§1. The first order equation 
1.1 Introduction. The linear differential equation of the first. 
order with constant coefficients is the simplest of the equations 
which will be treated in this book, and thorough familiarity 
with its properties will be absolutely essential in everything 
that follows. 
The equation is most conveniently written in the form 


Bt ay=f. (x) 


Here a is a constant (i.e. does not depend on x or y), fis a given 
function of x, and y is an unknown function of x which we have 
to find from (x) by ‘solving’ this equation. The equation is 
called ‘linear’ because it contains only a linear combination of 
y and dy/dx, as opposed to non-linear combinations such as 


of, and first order because it contains only the first deriva- 


tive dy/dx and not the higher derivatives d2y/dx?, . . . 
One instance of (1) will certainly be known to the reader: the 
equation 


dy _ 
a f- (2) 


Here we have to find a function y whose derivative equals the 
given function /, and this is simply the problem of integrating f. 
Thus y must be an indefinite integral of f: 


y=[seyarec, @) 


where C can be any constant. The general solution of (2) 
z 


LINEAR DIFFERENTIAL EQUATIONS 
therefore involves one arbitrary constant whose value is at our 
disposal. In practice it is often not the general solution which is 
needed, but rather some special solution, for instance that solu- 
tion for which y has some given value y, when x has some given 
value x». For example, if dy/dx=3x*—2 then y=x*—2x+C, 
but if we also require that y=1 when x=2 then we must take 
C=—3. Generally, if an initial condition: 
Y=Yo when x=%p, 

is imposed then the value of the arbitrary constant C in (3) can 
be determined. The required solution is then most simply 
written in terms of a definite integral of f, as 


I=Iot fi eae. «) 


Another familiar special case of (1) occurs when the right- 
hand member /(x) is zero, so that the equation reduces to 


B+ ay=o. (8) 


Dividing through by y and noting that 


xdy_d(logy) 
yan dx 
we can then write 
(log y)/dx=—a, 
whence 
log y=—ax-+C, y=C'e-™, 
where C is an arbitrary constant and C’ '. The preceding 
argument is, however, inadequate because it ignores the possi- 
bilities that y may be zero (then we must not divide by ¥) or, 
even worse, that y may be negative (then log y is not defined). 
It would not be hard to patch up the argument and to arrive 
at the complete solution to (5), but for us it will be better to 
make a fresh start: we shall give another method of treating (5) 
which not only leads to the complete solution of (5) but can 
also be used to solve the more general equation (1). 
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THE FIRST ORDER EQUATION 
1,2 The integrating factor. We now want to show that the 
general solution of (5) is 
y=Ke-™, (6) 


where K is an arbitrary constant. Now (6) is equivalent to 
ye"*=constant, 
which in turn is equivalent to 


d 
Z0e)=0, 


dy ney 
o( +4) =0; 
since ¢** is never zero, the last equation is equivalent to (5). 
Thus (5) and (6) are equivalent: every choice of K in (6) gives 
a solution of (5), and every solution of (5) has the form (6) for 
some value of K. 
The success of this method is due to the fact that 


dy 
a(t) 
is an ‘exact derivative’, i.e. is the derivative of ye. The same 
fact can be used to solve the general equation (1), 


Bs aymfa); 
if we multiply by e* we obtain 
Loe =erpt), 
and then we can integrate both sides to obtain 
yee [enfidec, 


moe enfide Com. (7) 
Formula (7) gives the general solution; if an initial condition 
has to be satisfied, say yy, when x=%, the formula would 
read 
gael ened tye ae. 6) 
% 
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LINEAR DIFFERENTIAL EQUATIONS 
In practice it is hardly worth while to try to memorize these 
two formulae: one need only remember the method by which 
they were obtained. 
Example 1: %_symsin * 
Multiply by e~* in order to bring the equation into the form 
Soe =e sin x, 
and then integrate to obtain 


ye = fo sin xdv+C 


a it Cc. 
= (cos * +3 sin x) +C. 
Hence the general solution is 
y=—= (cos x43 sin 2) +Ce™, 

where C is an arbitrary constant. 

Example 2: Solve Biyer for x>0, with the initial condition 
y=o when x=1. 
After multiplying by e* we get 

Joven axle, 
Integrating this and making use of the initial condition, we get 
2 
ye= J eetag 
1 
and hence the required solution is 
zy 
y = S-tetdg. 
x 

It will be noticed that our formula for the solution contains an 
integral which cannot be expressed in terms of ‘elementary’ func- 
tions (polynomials, exponentials, logarithms, sines and cosines). 


This will often happen and when it does we must be content to 
accept such formulae for a solution. In practical problems one is 
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THE FIRST ORDER EQUATION 
in any case likely to need numerical values for the solution, and 
from this standpoint the ‘known’ functions are those whose values 
have been tabulated: amongst these occur not only the ‘elemen- 
tary’ functions mentioned above but many others (including, as it 


z 
happens, the integral i) §-4e§d$ which occurs in Example 2). We 
1 


cannot of course expect that the solution will always involve 
integrals which have already been tabulated, and it may some- 
times be necessary to use numerical methods of integration (such 
as Simpson's Rule). 
The work of this section may conveniently be summarized 
in the following rule for solving the first order equation (x): 
Multiply by e* to bring the equation into the form 


d 
Se) =e7/). 


Then integrate both sides; initial conditions can be fitted 
at this stage by using appropriate definite integrals (as in 
Example 2). 
It is usual to call e* an integrating factor for equation (x) be- 
cause multiplication by this factor enables us to solve the 
equation by simple integration. 


1.3 The form of the general solution. We return briefly to 
the general formula (7), 


yao [enfildxs core, 
for the solution of dy/dx-+-ay=/(x). Taking C=o we see that 
goer |ere)ax 


is a particular solution, and that the general solution is obtained 
by adding Ce-**. Now Ce~® is the general solution of the 
reduced equation dy/dx-+-ay=0; it is usually called the comple- 
mentary function. Thus we can say: 
‘General solution=particular solution 
plus complementary function.’ 
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LINEAR DIFFERENTIAL EQUATIONS 
We shall find that a similar assertion may be made for linear 
equations of higher order and will then be very useful in finding 
the general solution. 


§ 2, The second order equation 
2.1. The reduced equation. The general second order equa- 
tion reads 
@y 
F954 Wy=fe), @) 


but we shall first treat the reduced equation in which the right- 
hand member is zero, i.e, 


ay & 
Bate tly=o. @ 
Here a@ and 6 are supposed to be constants; they must not 


depend on x. 
The first order equation 


+ ay=0, 


analogous to (2), has solutions Ce~**. We may therefore expect 

(2) to have exponential solutions also, and indeed we can easily 

find such solutions by substituting y=e** into (2). This gives 
Reh + aket*+-bet™=0, 

and so y=e"* will be a solution provided that 

k8-+-ak+-b=0. 

This quadratic has two roots, k, and ky, and so we obtain two 

solutions, e"* and e*#, of (2). Now whenever we know two 

solutions, say y, and ys, of (2), then Cyy,+4-Cyy, will also be a 

solution for arbitrary values of the constants C, and C,; for 

if y=Cyy4+Cyy_ then 


dy, dy 
wa ag? 
(coSot cet) +0( cca co2) Can tC) 
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THE SECOND ORDER EQUATION 
a’ ds 
-<( sd Pay) al? J a-+ty4) 
=C,(0)-+C,(0)=0. 
We have therefore found a solution 
y=Cyeh4-Cyeh 
which contains two arbitrary constants. This is quite satis- 
factory except when k, and k, happen to be equal; then we 
merely get y=(C,+-C,)e™*, i.e. an arbitrary multiple of e&* 
only, and we may expect that there is a second solution which 
we must still find. Also we cannot be sure (even when h, and kg 
are not equal) that we have found all possible solutions. We 
therefore make a fresh start, using a better method which will 
not only resolve our present difficulties but can also be used to 
deal with the general equation (1). 
The new method depends on ‘factorizing’ the left-hand side 

of (2) in a way which is perhaps made clearer ah tiie anew 


notation. Let us write ‘D’ for ‘d/dx’, so that wad, will be 
denoted by Dy, Then 

dty_d (dy 

Sg 2) ~PP»? 
and naturally we shall abbreviate D(Dy) to D*y. Equation (2) 
can now be written as 

(D*+-aD +-b)y=0, 

which at once suggests that we should factorize the expression 
in brackets as (D—h,)(D—h,), just as we would factorize 
k*+-ak-+-b as (k—h,)(k—k). Since ‘D’ is not a number but a 
symbol denoting ‘differentiation’, the reader will be quite 
justified in doubting whether such a factorization is legitimate. 
It is in fact legitimate because D satisfies the usual arithmetical 
laws relating to addition, subtraction and multiplication. To 
put it less abstractly: once we have agreed that D and D® shall 
stand for d/dx and d?/dx*, a reasonable interpretation for 
(D—hy)(D—hy)y is 


ree mm ae 
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dy ,y 
a Ryze t haley 


ay ay 
=D thi thy 


at, 

ar le ide 
because fy-}-ka=—a and hyky=b, hy and kg being the roots of 
k®-+-ak+b=0. 

‘We have now written (2) as (D—/,)(D—h,)y=o. Tf we put 
(D—h:)y=z, for the moment, this becomes (D—k,)2=0, or 
dz/dx—hyz=0. This gives z=Ce™, and now we have 

dy Be 
Y_ key=2=Cem*, 
zz! 
which is a first order linear equation for y. Solving it by the 
standard method we get 
Lae es 
shed 2) — Celie, 
wert 


ye*#=C j eta Wbedy LC, 


If ky-bhp, the integral on the Tight is 


© tts—w0e; 
Ry—ky 


if ky=hy, it is simply | Cdx—Cx. We have now proved that 
the most general solution of (2) is given by 


eh Cleh® if hy h, 


c 
doa h— 
y=Cxe+Cee* if hy=hy 
—hy) it arbitrary constant: we can give it 
he eae he Sie His suttebly choosing C. i Thus the 
general solution of (2) is an arbitrary linear combination of 
eh and e& if ky+-he, 
e&* and xe™* if ky=h. 
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THE SECOND ORDER EQUATION 
Incidentally we have found the missing second solution when 
hy=hy? it is xeh*, 

One further point needs discussion: what happens if the roots 
of k®+-ak-+-b=0 are complex? In practice a and b will certainly 
be real, so that 4, and , will be complex conjugates, say 

k=a+io, k,=a—iw. 
The corresponding solutions will be 
els tiole — ¢(cos wxi sin wx). 
Any linear combination of e+) and e(¢- can therefore be 
rewritten as a combination of e* cos wx and e* sin wx, say 
&(A cos wx-+-B sin wx); 
it is usually preferable to write the solution in this form because 
the constants A and B will then be real for a real-valued 
solution. 


Example 1: (D?42D+2)y=0. 
Here k*+2k +2 =0, so (k+1)?+1 =o andk+1= +i, ork=—144. 
Thus «= —1, w=1, and the solution is 
y=e*(A cos x+B sin x), 
Perhaps the reader who feels any doubt about the use made of 
complex numbers in this calculation may be slightly reassured 
when we have checked that this is a solution. We have 
Dy =—e-*(A cos x+B sin x) +-e~*(—A sin x+B cos #) 
=e-*[(B—A) cos x—(A +B) sin x), 
and similarly 
Dty=e*[—2B cos x +24 sin a}. 
Hence (D*+2D+2)y is a linear combination of e~*cos # and 
e~*sin x, with coefficients 
—2B +2(B—A) +24 =o, 
2A —2(A +B) +2B =o; 
that is to say, (D?+2D+2)y=0. 


2.2 The general equation. The ‘factorization’ method, used 
above to solve the reduced equation (D?+-aD-+-b)y=o, can be 
used equally well to solve the general equation (1): 


(D*+4aD+))y=/(x), 
E.D.E.—B 9 
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when the right-hand member /(x) is not zero. Putting (D—hy)y 
=2 as before, we now have 


(D—h,)z=flx). 
After solving this first order linear equation for z, we can then 
solve 
(D—ha)y= # 
for y. 


Unfortunately this method usually involves rather clumsy 
calculations because of the integrations needed for solving two 
first order linear equations, and the method is not recommended. 
for practical use. If we did use it, we should find that the solu- 
tion consists of two parts: one part, containing two arbitrary 
constants of integration, coincides with the solution of the 
reduced equation (see p. 8) and is called the complementary 
function; the other part depends on f(x) and is a particular 
solution. Thus the general solution of 

(D*+-aD-+b)y=f%) 
can be found by taking one particular solution and adding to it the 
rary function (i.e. the general solution of the reduced 
equation (D*+aD-+-b)y=0, which contains two arbitrary 
constants). 

In view of the practical importance of the above statement, 
we will give a direct verification. Let y, be a particular solution 
of (1) and let yy be any solution of the reduced equation. Then 
if y=9,-+y2 we have 

(D*+-aD +-b)y=(D*-+-aD-+0)y,+(D?-+aD-+0)y2 
=f(x)+0=/(%), 
so that y is also a solution of (1). Conversely, if y is any solution 
of (x) and if we put y—y,=Y2, then 
(D?-4-aD-+b)yy=(D*?-+-aD-+b)y—(D*+-aD+))y, 
=f) Sa) =0, 
so that y» satisfies the reduced equation; this shows that y is 
necessarily of the form +2. 

Now it is often possible to find a particular solution fairly 
easily without having to use the general method explained at 
the beginning of this section. In particular, this can always be 

ro 
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done when f(x) is a polynomial, an exponential, or a combination 
of sines and cosines. The next three sections will show how it 
can be done in these three special cases. 


2.3 Particular solution: polynomial f(x). If f(x) is a poly- 
nomial, we can always find a particular solution which is also 
a polynomial by substituting into the equation a trial poly- 
nomial (with undetermined coefficients) and then equating 
coefficients. 


Example 2: (D?-7D —5)y =4° —1. 
we try a particular solution of the form y=px*+g*-+-r4-+s. 
Then 
Dy =3px*+29%+1, Dy =6px +29, 
and we must have 
6px +-2q—7(3px? 4- 29% +7) —5(px* +9" +14 +5) =x8—1. 
Equating the coefficients of x, #*, * and 1, we find that 


I 
petted hte 


ar ar 
—21p —5q =0, = ——p=—} 
Msonim she 


6 204 324. 
149-57 =0, ———— = 57, r= — 
6p —149— 5 Sag ot tae 


mason a1, 2428 2603 
2q—7%—5s ag t x25 te SF eee 
Hence a particular solution is 

I= ghe(—125%" + 525%" —1620% +2603). 


Notice that the trial solution was taken to be a polynomial 
of the same degree as the right-hand side x*—z of the equation. 
‘This will suffice unless the coefficient b (of y on the left-hand 
side) is zero; in this case the degree of the trial polynomial must 
be increased by 1. The reason for this is that when 6=o the 
constant term s in the trial solution y= . . . +-s will disappear 
when we form (D?+-aD-+-o)y, and will not be available when 
we equate coefficients. We therefore need another coefficient 

Ir 
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at the start of the trial solution. If a, the coefficient of Dy, is 
also zero then the equation is simply D*y=/(x) and can be 
solved by integrating twice. 


2.4 Particular solution: exponential f(*). Another case in 
which the form of a particular solution can be predicted is 
when /(x) is an exponential, so that the equation reads 
(D?+-aD-+b)y=Mem, 
In general one can then find a particular solution of the form 
y=Ae™, and the constant A can be found by substituting into 
the equation. We must have 
Amen +a. Ame™ 4-5. Ae™™=Me™, 
A(m?+-am-+-b)=M. 
This fixes A, unless m?+-am-+-b=o. In this case, k=m is a root 
of the equation k?+-ak-++b=o: thus e™* occurs in the comple- 
mentary function, and naturally the result of substituting 
y=Ae™ into the left-hand side is zero. We must then try 
‘4y=Axe™ and find A as before by substituting. Finally it may 
happen that both e* and xe" already occur in the comple- 
mentary function, i.e. that k=m is a double root of k®-+-ak+-b 
=o. We must then try y=Ax*e™* and find A by substituting. 
(We have not attempted to explain why multiples of xe™* or 
x%em® should occur in the exceptional cases mentioned above. 
The reader may either accept this as an experimental fact or 
may derive it for himself by applying the general method of 
solution explained in § 2.2.) 
Example 3: (D?—3D +2)y =5e™. 
Substituting y=Ae™, we get 
oAe*—3.3Ae% +2Ae = 5e™, 


2Ae™ = 56%, A -2. 


Hence yaden is a particular solution, 


Example 4: (D*—1)y =3e~*. 
Here e~* already occurs in the complementary function; it is 
Iz 
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useless to substitute Ae~*, and we must try Axe~* instead. If 
y=Axe-®, then 
Dy=A(—x-+1)e~*, Dty=A(x—2)e*, 
So we must have 
A(x —2)e*—Axe~* = 367%, 
The terms in *e~* cancel as they should (this is always a useful 
check), and we are left with 


—2Ac*=3e-%, A= 3 
so that y= Sue is a particular solution, 


Example 5: (D*—2D +1)y =e. 
Here k*—2k +-1=0 has k=1 as a double root, and therefore both 
e® and xe® occur in the complementary function. Thus we have 
the worst possible case and must try y=Ax%e*. This gives 


Dy =A(x3--2x)e%, Dty =A(a#*+-4x +2)e%. 
Substituting and cancelling the common factor e*, we get 
A(x*+-4x-+2)—24 (x2+24) +A? =1, 


The terms in x* and # cancel (again a useful check), and this 
simply leaves 2A =1, A=}; y=}+*e* is a particular solution. 


2.5 Particular solution: trigonometric f(x). Finally, we 
deal with the case when /(x) is a multiple of cos nx or sin mx, 
or more generally of the form 

S(®)=P cos nx+-Q sin nx. 
There are several methods, of which we give two. First, we 
may try a solution of the same form as f(x), namely 

y=A cos nx+B sin nx, 

and find A and B by substituting in the equation. 


Example 6; (D*—2D—5)y =2 cos 34 — sin 3%, 
Substituting y=A cos 3x-+B sin 3x, we have 
(—94 cos 3¢—9B sin 3x) —2(—34 sin 3x +3B cos34) 
—5(A cos 3% 4B sin 3x) =2 cos 3x— sin 3%. 
13 
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Equate the coefficients of cos 3 and sin 3% to zero; this gives 


—144 —6B=2, 64 —14B =—1, 
and hence 


(142+6%)A =2(—14)—1(6), 2324 = —34, A=—17/116, 
(14 +-6*)B =2(—6)—1(—14), 232B=2, B=1/116. 
Thus we have as a particular solution 


T . 
Y= Fz™l —17 008 34+ sin 34). 


The preceding method will fail when cos mx and sin nx 
already occur in the complementary function, i.e. when the 
equation is 

(D?+-n*)y=P cos nx+Q sin nx. 
The trial solution must then taken to be 
y=x(A cos nx-+B sin nx); 
see Problem 12. 
A second method for finding a particular solution of 
(D?-+aD-+)y=P cos nx-+Q sin nx 
is based on the fact that cos mx and sin mx can be expressed in 
terms of complex exponentials. For instance cos nx and sin nx 
are the real and imaginary parts of e*, and more generally 
P cos nx+-Q sin nx is the real part of 


(P—iQ)(cos nx-+-i sin nx)=(P—iQ)e™*. 
(We shall assume that P and Q, as well as the coefficients a and 
b, are real.) Now suppose that we have found a complex-valued 


solution y of 

(D*+-aD-+b)y=f(x), 
where f(x)=/\(x)+i/2(x) is also complex-valued. Writing 
y=,-+ 472, where y, and y, are the real and imaginary parts of 
y, we have 


(D*-+-aD-+b)y,+i(D*+-aD +)y,=f,(2) +#fa(x)- 
Equating real parts, 
(D*+4aD+b)y,=/(x), 
and equating imaginary parts, ’ 
(D*-+-4aD+b)y.=fxl2)- 
14 
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Thus if y is a particular solution when the right-hand side is 
f(x), then Ry and Iy are solutions when the right-hand side is 
Rf(x) and If(x), respectively. To find a solution when the right- 
hand side is K cos mx or K sin nx, we may therefore first find 
a solution when the right-hand side is Ke, and then take its 
real or imaginary part; more generally when we have P cos nx 
-+@Q sin nx on the right, replace it by (P—iQ)e"*, find a solu- 
tion, and then take its real part. 

We illustrate the method for the same equation as in Ex- 
ample 6. 


Example 7: (D*—2D—5)y =2 cos 3x—sin 3%. 
The right-hand side is the real part of (2 +-7)e*. To find a particu- 
lar solution when the right-hand side is (2+i)e%*, substitute 
y=Ae™*, Then after cancelling e* we get 
A{(3i)*—2(34) —5]=2 +4, 
24h (2+4)(14—6i) _ 34+2t « 
RET. MAREGT TT age Re 

Taking the real part of the corresponding solution Ae**, we 
obtain the required solution 


gon{ —Te7 +(e sei aia »*} 


oars 4 — sin 34) 
=~ Fy (17 cos 3¢— sin 34). 


® A similar method applies in the exceptional case when 
cos nx and sin #x already occur in the complementary function. 
Example 8: (D? +-4)y =3 sin 2%. 
Here we shall first find a solution when the right-hand side is 
3e” and then take its imaginary part. Since e¥ already occurs 
in the complementary function we must take a trial solution of 
the form y=Axe"?, This gives 
Dy =A (2ix+1)e™*, Dty =A(—4x+4i)e¥, 

Substituting and cancelling e*, 


P Sie S> 
A(- 4=3; Anta —3i, 
(-44-+4)+444=3 ithstnrb 


Hence the required solution is 
y=1{—Jixe™} = — 4x c0s 2%. 
15 
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2.6 Particular solution: some further cases. The methods 
of the preceding three sections can be extended to deal with 
more complicated forms of right-hand side f(x). We shall 
merely state what procedures can be used: examples can be 
found amongst the problems at the end of this chapter. 
(i) If f(x) has the form 


e™*(P cos nx+Q sin nx), 
either substitute a trial solution of similar form, viz. 
y=e"*(A cos nx-+-B sin nx), 
or replace the right-hand side by 
(P—ignemtinw, 


find a solution, and take its real part. As usual a factor x must 
be inserted if e™* cos nx and e”* sin nx occur in the comple- 
mentary function. Unless the reader has strong objections to 
the use of complex numbers, he is advised to use the second of 
the above methods. 

(ii) If f(x) consists of a polynomial multiplied by an expon- 
ential factor, try a solution of the same form; the polynomial 
occurring in the trial solution should have the same degree as 
that occurring in f(x), but its degree must be increased in 
exceptional cases. For instance suppose that f(x)=(x?+-1)e™*. 
Then take y=(px?-+-¢x-+-7)e** as trial solution; but if e”* occurs 
in the complementary function take y=(px*-+-gx*+-rx)e™, and 
if xe* also occurs take y=(px4+-gx*+-1x")e™*. 

(iii) Similarly, if f(x) is a polynomial multiplied by (P cos nx 
+-Q sin nx), either try a solution of the same form or (better) 
replace the right-hand side by the same polynomial multiplied 
by (P—#Q)e'"*, find a solution as in (ii) above, and then take 
its real part. 

(iv) It may happen that f(x) is the sum of several terms, 
each of which could be treated by substituting an appropriate 
trial solution. Clearly we can deal with such an f(x) by sub- 
stituting an appropriate combination of terms. For instance if 
S(x)=e*—2e-%*-+x%, substitute y=Ae*+Be-*+-px*+-gx-+r 
and hence determine A, B, p, g andr. 

16 
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2.7 Arbitrary constants and initial conditions. The general 
solution of an equation 
(D?+-aD-+b)y=f(x) 

has the form y=¥,-+-4/, where , is a particular solution and yg 
is the complementary function containing ‘wo arbitrary con- 
stants, One often wants not the general solution, but rather the 
special solution which satisfies given initial conditions: y and 
Dy are to have given values at some given value of x, The 
appropriate values of the arbitrary constants can then always 
be found by substitution, as in the example below. 

Example 9: 

Find the solution of (D?—4D—5)y=% for which y=1 and 

Dy=—1 when x=o, 
The quadratic k?—4k—5 =o has roots k=5, —1,so that the com- 
plementary function is 4e%*+Be-*, As a particular solution, try 
y=px*+qx+r. This gives 

2p —4(apx +4) —5(ba* +92 +9) =a%; 


equating coefficients, 
—sP=1, —8p—59=0, 2b—4q—57=0, 
and hence 
pats q=— p= 5; sr—2p—q=-2 2, r= 
5 5 25 5 25 125 


So the general solution is 
y= 2994 by 42 4 sets Be-*, 
5 25 «125 
For this solution we have 
y= +448, 
é ? when #=0. 
Dy aces —B, 


To fit the initial conditions (y=1, Dy=—1, when x=o), we must 
have 


— +A+B=1, 
125 
rey | —B=-1, 
25 
17 
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Adding, we get +64 =0,4 == then 


Bartsdtnat epee 


75 27 3 
OT ars alae 
arts 4 
Bg eg ae —e¥+—e*, 
SE Syne ias wean a yse Sg 


The reader will notice that even a simple example may in- 
volve quite heavy calculations because one has first to find the 
general solution and then to fit the arbitrary constants to the 
given initial conditions. The arithmetical work can be lightened 
by using the operational method, to be presented in Ch. IT; this 
method is specially designed for finding the solution satisfying 
given initial conditions, and at the same time will also give the 
general solution if this is wanted. 


2.8 Recapitulation. We shall now give a brief summary of the 
methods developed in the preceding sections for finding the 
general solution of 
(D*4-aD-+))y=f). 
The reader will recall that the general solution is the sum of any 
one particular solution and of the complementary function; the 
latter is the solution of the ‘reduced equation’ 
(D?+-aD+b)y=o 
and contains two arbitrary constants. The following routine for 
finding the general solution is recommended: 
(i) Find the roots, ky and kg, of k*+-ak-+-b=0. 

(ii) Write down the complementary function ; this is 
Cye**-+-C,e' if ky and hy are real and unequal, 
Cye*+C,xe™* if k, and k, are real and equal, 

@*(C, cos wx-+-C, sin wx) if ky and ke are complex con- 
jugates aio; 
C, and C, are arbitrary constants. 
(iii) Find a particular solution. This can be done by substitut- 
18 
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ing a trial solution with undetermined coefficients in the 
Sollowing cases. 

(iii) f(x) is a polynomial of degree n. Try a polynomial of the 
same degree, y=px"-+-qx"+- ... ; if b=0, try a poly- 
nomial of degree (n-+-1), y=px"t14-9x"+ .... 

(iid) f(x) an exponential, say f(x)=Me. Try y=Ae™; this 
succeeds if m does not equal k, or ky and gives the partiou- 
lar solution 

M 
Tae 7 ie 
If ky=phg and m equals ky or kg, i.¢. if e™ occurs in the 
complementary function but xe"* does not, try y= Axe”, 
If ky=kg=m, i.e. if both e™* and xe™* occur in the com- 
plementary function, try y=Ax*e™. 

(iiic) f(x) is a@ combination of cos nx and sin nx. If f(x)= 
K cosnx (or K sinnx), replace f(x) by Ke, find a solution 
as in (iid) and then take its real part (or imaginary 
part). If f(x)=P cos nx+Q sin nx, replace f(x) by 
(P—iQ)e™, find a solution as in (jib) and then take its 
real part. (For another method, not using complex 
numbers, see § 2.5.) 

(iiid) f(x) is a polynomial multiplied by "= or by P cos nx 
+Q sin nx: see § 2.6. 

(iiie) If f(x) is not of one of the above types (nor a sum of several 
terms, each of one of these types) use the ‘factorization’ 
method, § 2.2: see Example 0 below. 

(iv) Add the aps tarsi fining found in step (ii), to the 
particular solution. If initial conditions have to be satis- 
fied, the arbitrary constants can now be determined, asin 
Example 9 (§ 2.7). 

The following example illustrates case (iiie) above, when a 
particular solution can not easily be found by guessing its form 
and substituting a trial solution. 

Example 10: Find a particular solution of 

(D?=—2D +1)y=1/(1+e*). 
19 
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Since the left-hand side can be written as (D—1)*y, we shall put 
(D—1)y =z, so that 
(D—1)z=1/(1 +e). 
Multiply this first order equation by the integrating factor e~*, 
which gives D(ze-*) =e~*/(1 -+e*) =e-**/(1-+-e~#), Hence 
Pad 
eae | * ae 
ze’ IF ae aie. 
Putting e~*=u, the integral becomes 
u I 
— | due — | (1 -—— )dus —u + 1 ; 
{ at ‘G : +) 1 log (1-41) 
we may ignore the constant of integration because we are only 
looking for a particular solution. So we may take 
ae-* = —e-* + log (1-+e-*), 
Now we must solve (D—1)y =z, or D(ye-*) =ze~#. 
This gives 


yom? me foe (r+e*)dx, 


yortet [log (1+0-*)dx, 


Note that this method will always require ¢wo integrations. 
In the above example, the first integration could be carried out 
explicitly (but the second could not); in general, however, 
neither integration can be done explicitly so that one obtains 
a rather awkward formula containing a repeated integral. 
Once again the use of the operational method (Ch, II) will 
improve matters because it gives solutions involving at worst 
single (and not repeated) integrals. 


§ 3. Equations of higher order and systems of first 
order equations 


3.1 The nt order equation. The general linear equation 
with constant coefficients, 
dry, dtty 
Pe ara at tay =f), (a) 
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can be treated by straightforward extensions of the methods 
used for second order equations. We again associate with (z) 
the reduced equation 


dty es 
mt +e +any=0, (2) 


and begin by treating this. It will be helpful to extend the ‘D’ 
notation by writing Dy for d*y/dx®, and so on. 

Our knowledge of second order equations suggests that we 
start by looking for exponential solutions of (2). We find at 
once, by substituting in (2), that y=e" is a solution provided 
that 


Raho}... + ay=0, (3) 
In general this equation has # distinct roots hy, ky,..., Rn» 
and then we get m distinct solutions e,... , e%*, When (3) 


has multiple roots, however, we get less than the expected 
quota of » distinct solutions and we must look for further solu- 
tions. Suppose then that k=/y is an r-fold root of (3); this 
means that the left-hand side of (3) can be written as 
P(k)(k—h,)", where P(#) is a polynomial which does not contain 
the factor k—/,. We now assert that the 7 functions 

eh, xebm, ete, |. | yf -Igkie (4) 
are all solutions of (2). To prove this we write (2) as 

(D"+-4,D""4- ... +-a,)y=0 
and then, by factorizing the left-hand side, write this as 

P(D)(D—hy)’'y=0; 
it will now suffice to prove that the r functions in (4) all satisfy 
the simpler equation (D—k,)"'y=o0. Now these functions all 
have the form e**#, where is a power of x, and we have 
(D—h,) (e*u) =D(e*u) —hye"u=eDu, 
(D—Ry)*(e#u) =(D—h,) (e"Du) =e%*D (Du) =e*D2u; 
by repeating this argument we finally obtain 
(D—hy)* (eu) =e%*D'u=0 
because D'u=o when u=1, x, x*,..., ‘1, This completes 
the proof that the functions (4) are solutions of (2). When we 
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have treated each root of (3) in this way, allowing for its multi- 
plicity, we can obtain the full quota of distinct solutions. 
Let the roots be hy, ka, . . . , 4, with multiplicities 7,, 7, .. . , 75, 
so that 7;+7,-+ ... +-7;=n. Then the » functions 

Oe, KOPP, IMIG, Oe, Stet; (5) 

sie 30, xef,..., #9 ee 

are solutions of (2). It follows that any linear combination of 
these functions is also a solution, and conversely it can be 
shown that every solution of (2) is a linear combination of the 
functions (5) (a proof of this last fact will be given later: see 
p. 63). Thus 


the general solution of the reduced equation (2) is an arbitrary 
linear combination of the n functions (5), where hy, .. . , ky are 
the distinct roots of (3), with multiplicities r,, .. . , Ty 


When some roots of (3) are complex, they will occur in pairs 
a-Liw (assuming that the coefficients a,, . . . ,4, are real), and 
one may then replace combinations of e+) and e*-* by 
combinations of ¢** cos wx and e** sin wx. 

Next it can be shown, by an t similar to that used 
on p. 10, that the general solution of the full equation (r) is the 
sum of any one particular solution and of the complementary 
function (general solution of the reduced equation), the latter 
involving » arbitrary constants. Particular solutions can be 
found whenever /(x) is a polynomial, exponential, cosine or 
sine, by substituting a trial solution of suitable form. The usual 
troubles occur when f(x)=Me™ and k=m is a root of equa- 
tion (3), ie. when e”* already occurs in the complementary 
function; if k=m is a root of multiplicity 7, one must take a 
trial solution of the form Ax‘e™*, 

It should hardly be necessary to give a detailed set of instruc- 
tions (like that on p. 18 for second order equations), but we will 
illustrate some typical points by an example. 


Example 1: Find the general solution of (D*—1)y=sin x. 
The roots of k?—1 =o are k= 4%, +1,so that the complementary 
function is an arbitrary linear combination of cos x, sin ¥, e* and 
e-*, To find a particular solution, replace the right-hand side by e* 
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(and take the imaginary part of the resulting particular solution). 
Since k =i is a root of kt—1=0, we must try y=A-e, This gives 
+-4i3)e8 = A(x —gi)e®, 
iy = —4ide®, 
and we must take —4id =1, A =ji. The particular solution is 
I{jixe} =jacosx, 

and the general solution is 

y =}scosx +C,cosx +C,sine +Cye?+Cye-*. 


3.2 First order systems. The methods of the present chapter 
can be adapted to problems involving simultaneous differential 
equations for several unknowns y, z,.. . , but the arithmetical 
work tends to become excessively clumsy. The operational 
method (Ch. I) is more efficient and systematic, and is there- 
fore to be preferred in all except the very simplest problems. 
Accordingly we content ourselves here with the simplest 
possible case: a pair of first order equations for two unknowns, 
y and z, 


dy =o, # ” 
Btw the=0, = +7y-+02=0, (6) 
with zero right-hand sides. 
We shall look for solutions of (6) in which y and z are both 
multiples of the sane exponential, say 
y=Ae™, 2—=Be', (7) 
Naturally we want a solution other than the trivial one y=z=o0 
so A and B should not both be zero. When we substitute into 
(6), we find after cancelling e** that 
(k+-0)A+BB=o, yA+(k-+6)B=0. (8) 
From (8) we obtain 
[(e-+-2) (+6) —By]A=[(k-+2)(k+6)—By]B=0 
so that A and B will both be zero unless f satisfies the quad- 


ratic equation 
(h-+-a)(k-+6)—fy=o. (9) 


When k does satisfy (9), then the left-hand members in (8) are 
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proportional and either equation in (8) determines the ratio 
A: B uniquely, (Those readers who are familiar with deter- 
minants will recognize (9), written in the form 

kta Bp 

h+6 

as a necessary and sufficient condition for equations (8) to have 
a non-trivial solution.) 

To avoid further complications we will now assume that the 
quadratic (9) has distinct roots, k, and k,. Taking k=, and 
k=h, in (8), ratios A, : By and Ag : By can be found such that 

y=Aye™*, z=Bye* and y=A ze", 2=Byeh* 
are solutions of (6). We can then form further solutions by 
taking arbitrary multiplies of these solutions and adding them 
together, and it can be shown that every solution of (6) can be 
obtained in this way. 

Example 2: 


=0, 


dy a NI} 
ae 3) F28=0, =F H4y—2=0. 


The quadratic (9) is (k—3)(k—1)—2:4=0, k*-4k—5=0, with 
roots k=5, —1, The equations (8) are 
(k—3)A -+2B <0, 44 +(k—1)B=o. 
When k=5, either equation gives A+B=o, and when k=—1 
either equation gives 2A —B =o. Thus we can take 
yue™, z= —e; 
or yre, s=2e-*; 
from these we can build up the general solution 
y=Cye™4-Cye®, z= —Cye®*--2C ye. 


3.3 Arbitrary constants and initial conditions. The general 
solution to the #'" order equation (1), p. 20, contains » arbi- 
trary constants. Their values can be determined if » suitably 
chosen conditions are imposed on the solution: usually these 
are initial conditions, and the » functions whose values are 
prescribed at some given point x=x, are y and its derivatives 
up to order »—1. We assert (without proof) that 
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the n™ order equation (1) has exactly one solution satisfying 
given initial conditions of the form y=yo, dy/dx=yy,..., 
arly /dxh-l=y, 1, when r=%q. 


This assertion will be proved in Ch, II, § 2.6 (though a differ- 
ent notation will there be used), 

Similarly the general solution of the pair of equations (6), 
Pp. 23, contains two arbitrary constants whose ‘Values can be 
determined when y and z are required to take given values 
Yo and % at x=%: 


the pair of first order equations (6) has exactly one solution 
satisfying given initial conditions 
=o, 2=%, When x=Xp, 


This is a very special case of a general fact about sets of first 
order equations which will also be established in Ch. II, § 2.6. 


PROBLEMS FOR CHAPTER I 


[y’ will denote dy/dx; Dy, D*y, . . . will denote dy/dx, 
d’y/dx*,... Unless pints stated, the general solution 
should be found.] 

1. y'—y=a?, 2. y' ye", y= 
3. y'+-2y=cos x, 4, y'+y 
5. (D®4-4D-+-5)y=2e7*, 6. (D8. . 
7. ras x-+sin x, 8. (D®—4)y=x*—3%—4. 
9. (D?4-2D-+1)y=xe~*, 10. (D°+-2D+5)y=x sin x. 
11. (D*—2D—3)y=0; y=2, Dy=—4, when x=0, 

12. Find a particular solution of MSc oat iti 
substituting the trial solution y=x(A cos x+-B sin x). 

13. eae oe 14, (D'—3D?+-D-+-5)y=100 cos 34. 

15. (D§—2D?+-1)y=x8. 

16. Dy=2y—z, Dz=3y—22; y=2, z=0, when x=0. 

17. Dy+y+z2=0, Da aese! 
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CHAPTER TWO 
The Operational Method 


§ 1. Preliminary discussion of the method 
1.1 The operator Q. The purpose of the operational method 
has already been indicated in Ch. I: it is to find solutions of 
linear differential equations with constant coefficients which 
satisfy specified initial conditions. 

We shall depart from our previous notations: the unknown 
will now be called x, and the independent variable ¢; the initial 
conditions will be imposed at ‘=o, This notation is in common 
use because it fits in with many applications in which ¢ is a time 
variable and x describes the state of some physical system 
whose initial state (at time 0) is given. 

To introduce the operational method we consider the initial 


value problem 
+ axef(), } 


4=%_ When t=0, 
for a first order equation. The solution is given by 


Sue) fter%, 
x(pe-e*—x(0)=[" Merde, 
° 


aang | ett-ofinae @) 


Thus (1) can be solved by an integration from 0 to ¢; but it can 

also be stated in a form involving only such integrations, with- 

out any differentiations. To see this, integrate each term in 
dx 

Tama) 
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from 0 to ¢, This gives 


x0) af ac)ar= far. (3) 


Conversely if x satisfies (3), then x(0)=x, because the integrals 
occurring in (3) are zero when t=0, and also by differentiating 
(3) we get 


Fall) —axl=/10. 


Thus the initial value problem (x) is equivalent to the ‘integral 
equation’ (3). The operational method takes (3) as its starting 
point, because this single equation incorporates the initial value 
%q Whilst in (x) the initial condition had to be stated separately. 

The operation ‘integrate from o to ¢’ will play an essential 
part from now on, and so we shall use a special notation for it. 
If f is any function of ¢, we shall denote by Qf that function 
whose value at ¢ is 


t 
on= | ‘see, @ 
and we shall speak of Q as the operator which produces Qf when 
applied to f. In this notation, (3) may be written as 
%—Xy—a0x=Of, 
(1—aQ)x=2%9+-O/, (5) 


where (1—aQ)x naturally stands for x—aQx. It is now tempt- 
ing to write down a solution to (5), namely 
x 


Q 
atta ©) 
regardless of the fact that 1—aQ is not a number but an opera- 


xm, 


tor so that the meaning of such expressions as 5 is not 


clear at present. However, we can tentatively identify the two 
terms on the right in (6) by comparing with the known solution 
(2); such a comparison suggests that we should interpret (6) 
by means of the rules 
I 


= a ” 
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(= J : et-Of(a) dr. 8) 


From Bip noe cea era 


, although this is not so evident 
because in (7) it omar only on the constant function x. 


operator; so ronst = 


1,2 Formal calculations with Q. So far no advantage has 
been gained from the introduction of Q, since we already knew 
the solution to problem (r). Let us now look at an initial value 
problem for a second order equation: 


dx dx 
mtg tS, ©) 


x=%p, Han, when t=0- 


We convert this into a form involving integrations only by 
integrating twice. The first integration gives 


% x, 4alx—x) +00x=0f, 


Sp artiQras, tay tO ; (x0) 
the second integration gives 
#—Ho-+-aQx-+b0(Qx)=O(x, +4%) + 0(Q/), 
x+aQx+bQ*x=%)+a0%+On+O%F, (x1) 


where we have written Q?f for Q(Qf) so that Q* denotes re- 
peated integration. Now suppose conversely that x satisfies 
(tz). Then «=x, when t=o0, because all terms in (11) which 
involve Q or 9? are zero when ‘=o. Next, by differentiating (11) 
we get back to (ro), and on putting ‘=o in (10) we find that 
Bax, when ¢=o. Finally, by differentiating (10) we get back 
to (9). Thus we have again converted the initial value problem 
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(9) into an equivalent single equation, (11), which im 
integrations only and incorporates the initial values x9 
We now write the left-hand side of (11) as 


(t4+-4Q+00")x 


(z—20)(z—nQ)x, 
where 4 and yv are the roots of 4*-++-ak-++-b=0; this brings (11) 
to the form 
(2—4Q)(1—nQ)x=%q + aQ% + Om + OF. 

Equations of this kind can be solved quite easily, provided that 
we are prepared to manipulate expressions involving Q rather 
freely; the following example illustrates what kind of mani- 
pulations are likely to occur. 


Example 1: 


and factorize this as 


dt dx 
= —3G tefl), 


=I, oo, when t=o. 


By integrating twice we get (as in (11)) the equivalent equation 
(1-32 +20*)a =1 —30(1) +O(—1) +0%, 
(1-Q)(1 -20)* = (1-49) +0. 
Proceeding quite formally, we get 
Parent Aa” a ne; 
(=a 20)" * (=O) =a) 
To: Geal, wi the dia, Seees 0 the right, we note that we can 


evaluate seri and = a by using formula (7), p. 27. We shall 


therefore evaluate 


G ae ne by putting it into partial frac- 


tions; thus 
1-40 3 2 
90-29)" = ig)! atic 
a 
We can deal similarly with ——y 9) 
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| EUR ye plan a a 
6-0-2) € apt) 9! tah 
the last expression can now be evaluated from formula (8), p. 28. 
So we arrive at the tentative solution 


eaget-f2ety f “(etn —elt-myf(eiar, 
0 


Some readers may find it interesting to verify that this is a 
solution; to do this, they may prefer to write the last term on 
the right as 


t t 
of e-*f(r)dr—et e-ef(a)ar. 
0 0 


Naturally we must still show that the treatment to which we 
subjected Q in the above example is legitimate. This will be 
done in the next few sections. However, some readers may be 
eager to learn the technique of the operational method as 
quickly as possible, and may be prepared to accept formal cal- 
culations with Q at face value. Such readers may, if they wish, 
pass straight on to § 1.6 in which the technique will be further 
developed. 


1,3 Operators. Whenever we have a rule for producing from any 
given function f a new function Af, we shall call A an operator. 
For example, the rules 


ar= foto. 


° 
Df)=f', 
define the operators Q (integration) and D (differentiation). By 
‘function’ we shall always mean a function of ¢ which possesses 
derivatives of every order, i.e. which can be differentiated as 
many times as we please. 
We call A a lineay operator if 


AGf+ug) =14f-+uAg, 
for any functions f, g and constants 4, 4; here A(4f+-g) stands 
for the result of applying the operator A to the function 4f+ug. 
Clearly Q and D, as defined above, are linear. 
If A and B are linear operators, they may be combined in two 
ways. First, we may take two constants « and #, and form the 
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linear combination «A +B. The effect of this operator on any 
function f is defined by 

(2A +BB)f=aAf+fBf. 
Secondly we may form the product AB. The effect of this operator 
on f is defined by 
AB{=A(Bf). 
Thus: 


to find ABf, first apply B to f, and then apply A to the vesult Bf. 


It is almost obvious that «4 +fB and AB are linear operators 
when A and B are linear, and we shall omit the formal verifica- 
tion of this fact. 

The definition of product will be illustrated by calculating QD 
and DQ. To find QDf, first apply D to f (i.e. differentiate) and then 
apply Q (i.e. integrate from o to #); 


anpi)= ['Ftede=s9 -S(0). (02) 

Similarly oe 
Dorin=F I Seidens 0. (03) 

0 


Notice that QD and DQ are not identical, and we must therefore 
take care to attend to the order of the factors when dealing with 
products of operators. When AB and BA are equal, we shall say 
that A and B commute. 

When one performs numerical calculations, one uses the laws 
of algebra quite freely. We will not give an exhaustive list of these 
Jaws, but recall several typical ones: 

a.b=b.a (‘commutative law of multiplication’), 
a.(b-+c) =a.b+a.c (‘distributive law’), 

a.(b.c) =(a.b).c (‘associative law of multiplication’), 
Analogous laws also hold for calculations with operators, with one 
important exception (already noted): the commutative law 
AB=BA does not always hold. All other laws relating to addition 
and multiplication continue to hold; for instance the distributive 
law, A(B +C) =AB-+AC, and the associative law, A(BC) =(AB)G, 
where A, B and C are linear operators. The associative law means 
that we can simply write ABC, without any brackets, for the 
operator whose effect on f is found by first applying C to f, then 

3r 


"a 


THE OPERATIONAL METHOD 
applying B to the result Cf, and finally applying A to the result 
Bf). 

In particular we may form the powers of a single operator A: 

A%=AA, A*=AAA,..., 
so that A" stands for ‘A applied m times in succession’. We can 
then form polynomials in A: 
CotGyA +0,A*+ ... +o_A", 
where ¢, denotes the operation of multiplying by cy. Any two such 
polynomials commute, so that polynomials in A may be multi- 
plied together or factorized just like ordinary polynomials, the 
order of the factors being immaterial. For instance 
(14-A)(1-2A) =(1-2A)(1-+A) =1-A —2A*, 


1.4 The inverse of an operator. Now consider an equation of 
the form 
Ax=g, 

where A is a linear operator, g is a given function and x is an 
unknown function. (Equations of this kind have already occurred 
several times, for instance (5) on p. 27 and (11) on p. 28.) Such 
an equation may have no solutions at all, or on the other hand 
may have infinitely many solutions. For example Qx=g, that is 


[ixereeat0, 
i 


has no solution if g(o) +0 because the left-hand side is necessarily 
zero when ¢=0; but Dx =g, that is 


dx 
Pricdod 
always has infinitely many solutions: 


t 
x()= g(t)dr+C. 
0 


However, suppose that Ax=g has exactly one solution for each 
given g. This solution will of course depend on g; we shall denote it 
by A~'g, and call the operator A-} the inverse of A. The operator 
A-1 is again linear, i.e. we have 

A614 162) =24-%9, + HAE 
To prove this, write x, =A~1g, and x,=A~1g,, so that Ax, =g, and 
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Ax,=g,. Then 
A (tay + rq) =1Ax, +HA%=]g, +H8y 
and therefore =x, +x, is a solution of the equation 
Ax=1g,+H8s 
But this equation has only one solution, namely 
#=A-'(0g, +082), 
and so 2%, +.*, must coincide with the above expression; in other 


words we have 
4A", + uA~'gg= AAG +82), 
which proves our assertion that A~? is linear, Next, we have 
A(A~'g) =8; 
this merely states the fact that x=A~g is a solution of Ar=g. 
But we also have 
A~"(Ag) =8, 


because = A~1(Ag) is by definition the solution of Ax = Ag; now 
x=g is clearly a solution of this equation and since the equation 
has only one solution, g must coincide with A~1(4g). We have 
now evaluated the product of A and A™, in either order, and our 
results can be written as 
AA-\=1, A-1A =1, 
where 1 denotes the ‘identity operator’ defined by 1g=g. Con- 
versely, suppose that we can find a linear operator B such that 
AB=1, BA=1, 
Then B is the inverse of 4; that is, the equation Ax =g has exactly 
one solution, given by +=Bg. To see this, note that x=Bg is a 
solution because 
A(Bg) =ABg=1g=g, 
and that it is the only solution because if Ay—=g then 
x =14=BAx=B(Aa) =Bg. 

To sum up the preceding discussion: 

If the equation Ax =g has exactly one solution, x=A-g, for each 

given g, then the operator A-* is linear, and AA~1=A—1A =1; 

‘AW will be called the inverse of A. 

Conversely if there exists a linear operator B such that AB =BA =1, 

then A has an inverse, namely A-?=B. 


To illustrate the definition of the inverse operator, let us first 
observe that neither Q nor D possesses an inverse; the equation 
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Qx=g sometimes has no solution, whilst Dx=g always has in- 
finitely many solutions. Next, we will show that 1—aQ (where @ 
is any constant) does have an inverse; this important fact forms 
the basis of the operational method. To establish it, we must look 
at the equation 
(t-aQ)x=g, 


that is . 
x(t -af a(t}dr=e(). 
0 


This is equivalent to 

# (t) —ax(t) =8"(t), *(0) =g(0), 
and we know from Ch. I (§ 1.2) that this initial value problem has 
exactly one solution: 


a(l)=e(o)ert+ ferreeee (4) 
0 


Thus (t—aQ)- exists, and x=(r —aQ)-%g can be calculated from 
formula (14). (At this point, we should recall the agreement made 
on p. 30 that all functions are to have derivatives of every 
order. Now 2(t), given by (14), is such a function; its first deriva- 
tive exists and equals ax(t) +g'(¢), its second derivative therefore 
exists and equals ax’(t)-+g’"(t), and so on.) 

In particular if we take g(!)=1, identically, we have g(o)=1 
and g’(z) =o; if we take g=Qf, then g(o) =o and g’(r) =f(z); from 
(t4) we therefore obtain the formulae 

(t-aQ)1 =e", (15) 


(—aay-40r0 = f ‘eot—0Y (har, (36) 
0 


These are precisely the formulae (7) and (8) which me tentatively 
wrote down in § 1.1, except that there we wrote ——— 


ana 2 
r= oe —aQ 
instead of (1—aQ)-! and (1—aQ)-9. 

1.5 Inverse of a product. When we come to deal with second 
order equations, as in § 1.2, we are led to equations of the form 
(1 +4Q+60)4 = (+ 6,0)1 + 0%. 

It will then be useful to know the inverse of 1 +aQ+bQ?; we shall 
see presently that this operator has an inverse because it can be 
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factorized as (1—4Q)(1—Q) and we already know the inverses 
of the factors 1—JQ and 1—Q. 

Consider then the following more general question; given that 
two operators A, B have inverses, does AB also have an inverse? 
To answer this question, we look at the equation 

ABs =g. 

This can be solved in two steps; first write the equation as 

A(Bx)=g, 
and solve for Bx. This gives 

Bz=A-'g, 
and now we can solve for x, obtaining 

x=B-}(A-1g) =B-1A-g, 
This is the only solution, and therefore AB does have an inverse, 
given by 
(4B)-2=B-14-1, 

We can also verify this by showing that the product of AB and 

B-1A-1, in either order, is 1. Thus 
(AB)(B.A-~) =A(BB-) A =A(1)A-=AA-51, 
(B-14-1)(AB) =B-(4-14)B =B-1(1)B =B—1B =1. 

Now suppose in addition that A and B commute. Then 

(BA)-*=A—B-1, but BA =AB so that we can now write 
(AB)-? =A-1B-1 =B-14-1, 
Thus: 
if A and B commute and have inverses A~! and B-}, then AB 
has inverse B-\A~1=A-1B-1, and so A-* and B-) also 
commute, 
In particular, if A has an inverse then A*=AA has inverse 
A-}A-1=(A-1)2, AS=A*A has inverse (4%)~14-1=(A-3)44-1 
= (A-*)', and generally 
(4")-13 (4-1), 
We shall therefore denote the inverse of A" by A~", and we then 
have A~"A"=A"A-" =1, Indeed we have now defined A’ for 
positive and negative integers r; if we define A®=1, then the ‘index 
law’ ATA*=A’** will hold for ally and s (positive, negative or zero). 

The following fact will also be useful: 

1 {if A has an inverse, and if A commutes with B, then A* also 
with B; that is, A-1B =BA-", 
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To prove this, take the equation AB=BA and multiply by A-* 
at the beginning and end: 
A-\(AB)A-1=A-\(BA)A~1, 
(4-14) BA-1 = A“1B(A A>»); 
since A~1A =A A~1=1, the last equation reduces to BA—!=A-1B, 
as required. 

The preceding work may seem remote from the practical prob- 
lem of solving differential equations, but in fact it will be very 
useful to us because we shall be dealing with products of operators 
in which the factors are all of the form 1—Q, 1—Q, .. , and so 
commute with each other, and also commute with Q. We shall 
therefore be able to use rules I and II above quite often, and by 
doing this we shall be able to handle calculations with Q just like 
ordinary (numerical) calculations. For example, II shows that 

(1-aQ)-19 =Q(1 —aQ)-, 
and I shows that if 1-+aQ+0Q=(1—AQ)(1—Q), the inverse of 
1+aQ-+bQ? exists and is given by 

(1 -+4Q +69*)—* = (1 —2Q)-(1 —4Q)-* = (1 —xQ)-*(1 —2Q)-2, 
Again, if we have an expression like 

(2-Q)-"(x —39)-"(t—Q)-"(t — 3) 
we can ‘cancel the common factor 1—3Q’ by using I to reverse the 
order of the two middle factors, which gives 


(t= Q)-(r-Q)-"(1 —3Q)-2(t —3Q) = ((t-Q)-2)*=(1—Q)-*. 
We can now afford to simplify our notation, without any risk 
of ambiguity, by writing 


for (t—aQ)-}, 
a a 
ao ~aQ)-19 =0(t-a0)-*, 


and so on; more generally if we have any expression 


$(Q) 


(1a, 
where 4(Q) is a polynomial in Q, it shall have the interpretation 
(1-a,Q)-". ‘j (1 4nQ)*9(Q) 
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and we may then alter the order of the factors in any way we 
please. If it so happens that ¢(Q) contains say (1 —a,Q) as a factor, 
we can cancel this factor against (1—a,Q)- by a suitable re- 
ordering which brings (1 —a,Q)-! and (r—a,Q) next to each other. 

It is now easy to justify ‘partial fraction expansions’ such as 
those which we used tentatively in Example 1 (p. 29). Thus, to 


—42_ 5 : 
put fe a (20) into partial fractions, write the numerator as 
1~4Q =3(1—20) —2(1—9); 
then. 
1-40 3(t—29) 2(t—Q) 3 2 


(1=Q)(x=20)" (r-Q)(1—20) (t—-O)(X—20) 1-9 1 =29" 
on cancelling the appropriate common factors. Similar justifica- 
tions can be given in more complicated examples; thus any 
rational function of Q, of the form 
Got OQ+ . + + ton"? 
« =a) - + (1=4nQ) 

may be put into partial fractions of the form 

A As An 

1a, *1=a,Q* +++ +700 

just as if Q were an ordinary (numerical) variable. (Of course the 
form of the expansion must be modified when the denominator 
(z—@,Q) . . . ({—a,Q) contains repeated factors; this point will be 
dealt with in the next section.) 


1.6 Partial fractions for inverses. We have now reached a 
point where we can justify any of the formal calculations with 
Q which we shall perform, and it is time to develop the tech- 
nique of such ere somewhat further. We already know 


how to evaluate — ay and > , but we shall also need 
the more general eeiont 
mei" grt 
G=egF® Gaagenl HOE AS) 
which we shall now evaluate. 
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First, using (7) and (8) from pp. 27-8, we have 


2 eas 6) I Q 
ana aalmag) a 


a f * alt—verd 
s t 
=e rdr=tett, 
o 
Clearly this argument can be repeated: 
e@ Q Q 
Ga (= rag 


a f } ot-ratdrmet avmte; 
0 0 2 


een (md ae male") 


tz2 
aa Caran 
° 
Proceeding in this way, we get the general result 
Leer 
(age al 7) 
for n=1, 2, 3,... ; it also holds for »=o when we make the 


usual definitions that Q°=1, =1 and o!=r. 
Notice that we have chosen to evaluate atone rather 


(t—aQ)* 
aor The latter is given by 


ra 


=e" alet=(r-+-at)e™ 
and is thus less simple than a Cis 
The second formula to be established reads 
nk 
mage= -f. Moya. (a8) 
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This holds when »=o, and we shall prove it by induction for 
n=I, 2, 3,.... Suppose that (18) has been established for =k. 


Then 
gus gi iy) ge 
G=agre! ~ a9 cane} 


gat = [ato Oeeatebae 


= ferroff" Ma —_ Preeayuyauhae. 


u=T 


u 


Tet 
In this repeated integral we can invert the order of the integra- 
tions, provided that we attend carefully to the limits of integra- 
tion, We have o<u<r<t, so if we do the t-integration first (keep- 
ing « fixed), then t varies from % to #; in the second integration, « 
varies from o to ¢, (This can also be seen by regarding the integral 
as a double integral, taken over the shaded triangle in the (r, 1)- 
plane in the diagram. If the w-integration is done first, goes 
from o to 7; if the t-integration is done first, r goes from % to é; 
in either case the second integration is from o to ¢.) Thus we get 


ff, catt-nfF, de eee 
Dees (r—u)* 
“fi: ” woffa 


« saewp, tA 
- | 0 de 
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This is (18) for n=k+1, except that the integration variable is 
called % instead of r, and the induction has been completed. 


Armed with formulae (17) and (18), we can now tackle any 


equation like 
(1+aQ+9*)x=(4,+0,0)1+0%, (19) 
and also similar but more complicated equations which arise 
from differential equations of higher than second order. Let 
(r+-4Q+60*)=(xt—4Q)(t—Q) ; then the solution to (19) is 
G+e20 oe 
= aaa , 
(= A9\0— nO)" = 20) = OF” 
and this can be evaluated by putting the right-hand side into 
suitable partial fractions. Example 1 (p. 29) shows how this is 
done when 4 and y are unequal; the following example shows 
how one proceeds when 4 and yu are equal. 
Example 2: 
a= Q* 
Evaluate 7=29+9"" Ks ago 
Here cig te —Q)*. We could express ——~ cS 


T 
ce es ap™ SER ul ae Vik oe an =O? 
rather than ———. tor ae So we write 2—Q as 2(1—Q)+Q, and then 
2-2 _ 21-9) +2 2 Q 
G=9F ~ @=9)* *W3=9" t=O 
=2e!+tet, by (17). 

The second term can be evaluated at once from (18), with n=1, 
and we therefore le 
on rte pelos l= (2 snes f' (t—z)el-*f(r)dr. 


Similar ven apply to the evaluation of the more general 


we have 


te aa ee 
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where $(Q) is a polynomial of degree less than n. Factorize 
1+4,Q+ ...+4,Q" as (I—/,Q)™... (I—A,Q)™2, and then 
express 

___ 4) 
(t-A,Q)™. . . (@—A,Q)™ 
as a combination of 
Q oe 


12,0’ (t—4,0)" ***’ (x—2,0)™ 
and similar termsinvolving 1—A,Q, . . . , I—A,Q; likewise write 
te Csi as a combination of 
I+... +a,Q" 
Q Qu 


1—A,Q' "*"’ (x—4,Q)™ 
and similar terms. The expressions (20) can then be evaluated 
by using formulae (17) and (18). 

In general the use of formula (18) will lead to integrals which 
cannot be further simplified. However, in practice the function 
f (which stands on the right-hand side of the differential equa- 
tion) is often a simple exponential, or polynomial, or a product 
of a polynomial and an exponential. The solution can then be 
written in a form which contains no integrals at all; the follow- 
ing example shows how and why this can be done. 


Example 3: 
Find the solution So al 
Sor which x=2 and So when t=o. 
The equivalent equation involving Q is found as on p. 28; it is 
(t-Q—29*)*=(2—20)1 +9%(3e~). 
‘We now use (17) to replace 3e~t by aot this gives 
(1-9-2992= (2-291 + 32, 
(r+Q)(r -2Q)x=2 = rad 
peek ae +0" 
GF O0—20)" 
E.D.E.—D 4r 
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ae 
(1+9)*(t—29) 


ides 
Ty oF 


Now write as a combination of —— 


Bi ey 
1-29" 


i ewe = 

1+Q@ (1+Q)? 1-20 
and finally evaluate this by using (17). The solution is 
#=e*—tet 0%, 


The above procedure can clearly be applied whenever /(¢) is 
a multiple of #e* or a linear combination of such multiples. 
Only formula (17) is needed, and the result will be of the same 
general form as f(t) (but with some extra exponential terms 
which belong to the complementary function). 


§ 2. Practical instructions for using the method 


2.1 The symbol p. The general principles underlying the 
operational method have now been established. For practical 
purposes, it is convenient to make a change of notation which 
will help to simplify the algebraic work involved in handling 
rational functions of Q. We introduce a new symbol, », make 
the substitution Q=/-", and then perform all our calculations 
with /, treating expressions in # according to the usual alge- 
braic laws. 


Inthis So ~ 


™ 
(r— Apne 
Thus the basic formula (17) from 


— Sana by ~ 


which simplifies to 
p- 38 becomes 


gt Aynee 
ym 
@- ia Sat 7) 
and similarly (18) becomes 
t (pm 
ganal= | Sow oeae. @) 
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The special cases 
Fa 4 
p A 
3 shone oe at I=te"; (3) 
Fn0= | feat, f= ‘ M—Yic)ae 
should be particularly Saige 
The typical expression 
te <= ey ey 
THa,0+ - tag" It... Hay 
now becomes 
Cob" IOP... +a, x 
a P+ Mpr +... +n ee iota 
ce) 
~Pagy ag 


The advantage of the latter form is that we can put C()/A(p) 
and 1/A() into partial fractions of the usual kind, containing 
terms in 

ee peti 

p—i (b—A)* ""** (6—2)™ 
for each factor (p—/)™ of A(p), and then 


C() 
Pry ta ap" 
will contain terms in 
p 
Bah Goat pale GwO 


which can be evaluated from formulae (1) and (2). Standard 
methods for finding partial fractions are therefore available; if 
we were to use Q we should have to modify these methods 
because we should need fractions of the less convenient form 


x Q a 
1—AQ' (x—10)" **** (@—10)™ 
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It should be stressed that the symbol p merely serves to 
simplify calculations, and that it is notan operator; the operator 
on which the method depends is p-!=Q. (If p were an opera- 
tor, it ought to be the inverse of Q; but, as we saw on p, 33, 
Q does not have an inverse.) The justification for using # lies 
in the fact that any valid calculation with Q can be ‘translated’ 
into an equivalent calculation with p, and vice versa; for 
instance the -calculation 


2p+X hee p 

‘a ee oe 
is equivalent to the ae Sie 

240 _Ht+Q)+H—)_ 1 4 


1-@ 1—@? 77-0" 
and we already know that such calculations with Q are per- 
mitted. 


2.2 Procedure for solving n‘® order equations. To deal 
with the initial value problem 
(D"+a,D"4+ ... +a)2=f(é), () 
x=%q, Dix, ..., D*-x=%__, when t=0, 
we integrate m times. This leads to an equivalent equation, in 
terms of Q, into which the initial values x», ..., %,_, have 
been incorporated, If we then put Q=p~? and multiply by p" 
in order to remove negative powers of /, we shall arrive at the 
following equation: 
(Prbaypr t+... +an)x 
mpl) LTP +E «by af) bm ("* 
TO ai oe beep ttm apfe © 
This will be called the operational. form of (4). The details of the 
calculation which leads to it will be given later (p. 62); we can 
check it for n=2 by putting Q=p- in (r2), p. 28, and multiply- 
ing by f°. (Note that @ and 4 in (rr) are now called a, and a,.) 
At present we want to show how (5) can most easily be memor- 
ized, and we recommend the following rule for finding the 
operational form of (4): 
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On the left-hand side of the differential equation, replace D by 

the symbol p; 

on the right-hand side, add to f(t) the polynomial in p which is 

obtained from 

rapt... be) (sot. + F8et)x 

by multiplying out and retaining only the positive powers of p. 
It should be noted that the polynomial on the right of (5) 
always contains / as a factor, and also that none of its terms 
contain the coefficient a,. 

The function /(¢) is often the sum of one or more terms of the 
type Afe*. Using (1), we can then write /(¢) as a sum of terms 


Ans 1; the resulting expression, when simplified, will 


ya’ 
have the form F(f)r and will be called the operational form of 
4). For example, #24 (3t—2)e~* has the SAD form 
2 p p  P(I— 2p), 
ga pe a Gat)" 
When we replace f(¢) by F(p)x in equation (5), we obtain 
A(p)x=(6C(6)+F (P))t- 
If we solve for x, writing the solution in the form 


Gp) 
*=p—1, 

H(t) 
we can then put the rational function G(p)/H(p) into partial 
fractions and evaluate x by using formula (r). It is important 
to remember that a factor p must be extracted before finding 
the partial fractions, pray we want to express x as a com- 
I 
bination of terms —~—x ( and not —— 
(e- = ( (p—A)" 
accounted for in the following rule for solving (4) when the 

operational form of f(t) is known: 

Write down the operational form (5) of (4), and at the same time 
replace f(t) by its operational form F(p)1. Solve for 3 obtaining 
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itin the form —=R(p)t where R(p) is a rational function. Put 
R(p) into partial fractions, and then evaluate x=pR(p)x from 
Sormula (x). 
Let us now illustrate the procedure by two examples. 
Example 1: Solve (D* +-2D* ~D —2)x =3e*, 
Sor the initial conditions xy=2, x;= —3, ¥,=1- 


The operational form of 3e% is 3-1; the positive powers of pin 
I 


wr+6t—o(2—3+ ) are 2p? +-p"—7p. Hence the operational 


Pp 
form of the problem is 
(et +29*—p—abe= (26*49*— np eat 
=7t30* 90" +170 
=e 
Factorizing p*+ap*—p—2 as (p—1)(p-+1)(P-+2), we get 
Haba —op417 
> @-Ne+e+2)—2) 
(545 gira 3/4 )s 


Pi ptt pa p=2 
(We have suppressed the calculations required for the partial 
fraction expansion; the technique of such calculations will be 
discussed in the next section.) Thus we have 
Pb 


7( 2 cD ft a (ale : 
Ute) AG ae 


PP ac 8 a LP 
hie ie we 


The reader will observe that the first three terms belong to the 
complementary function, and that the last term is a particular 
solution; he is also invited to check that the initial conditions are 
satisfied. 


Example 2: Solve (D?—7D —6)x=t*, 


Sor the initial conditions x,=x,=0. 
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Since 1# =r, the operational form of the problem is simply 
(6-76 +6)9— 
no other terms appear on the right because the initial values 
# and #, are zero. We now have 
ee Leak See 
Pp Pp P—1)(P—6) 
= (=2/5 1/549, 43/108 7/28 1/3), 
ee, Spero eat 


et dae 43,72 Tr 
o-(Eta teste ia): 


ee 
708 T1786" 

We now give a rule for solving (4) when the operational form of 
Sb is NOT known (i.e. when f(t) is not the sum of one or more 
terms Affe**): 


Write down the operational form 
A(p)x=pC (PI +f) 
of the problem. Then 
CP) 3 
tay apy 
This is to be evaluated by putting C(p)/A(p) and 1/A(p) into 
partial fractions, and then using formulae (1) and (2). 
Note that once again a factor p should be extracted before the 
first partial fraction expansion (but not before the second one). 
Example 3: Solve (D?—3D+2)x=f(0), 
for the initial conditions x)= —1, %,=5, ¥,=2. 
The operational form, found from the rule on p. 45 or from (5), is 


(p?—3P +2)*=(—1(6*—39) +5p*+20)t +h) 
=(—p*+5P*+5p)t +4). 
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Since p*—3p +2 =(p —1)*(p +-2), we obtain 
ies el PS SR, 
Pept) =p 


3 I —1/9 1/3 1/9 
“(pip as) (Gr pote 


t 
mztet—o-tt 4. [J —Zett- 4 2(¢—nyett—d 4 26-mtt-) | pea)ar. 
ol 9 3 9 


2.3 Some remarks on partial fractions. In the preceding 
examples we have deliberately suppressed the arithmetical 
work involved in finding the various partial fraction expan- 
sions. We must now warn the reader that this work usually 
forms a substantial part of the operational method of solution. 
It is therefore desirable to have a systematic technique for 
finding partial fractions; we shall now describe such a tech- 
nique, and shall use it in all future examples. 

Suppose then that G(s)/H(p) is a rational function, the 
degree of the numerator G(#) being less than the degree of the 
denominator H(p). Let p—A be a typical factor of H(p), with 
multiplicity m; thus H(p)=(p—2)™K(p) where K(f) does not 
contain the factor —A. Then there is an expansion 

G()__ An 

He) @—A* 
where the last dots represent similar terms corresponding to 
other factors of H(p). (We assume here that G(p)/H(p) can be 
expressed in the form (6); this fact is proved in textbooks on 
algebra, e.g. in Durell and Robson, Advanced Algebra, Vol. I1.) 
We want to find the coefficients A,,,..., A, in (6), and in 
doing this we are permitted to treat # as a numerical variable. 

When p—A is a simple factor of H(p) (i.e. when m=1) the 
coefficient A in the single corresponding term A/(/—A) can be 
written down at sight. For then, if we multiply the identity (6) 
by H(), we obtain 


G(p)=AK(p)+(P—A)L (0), 


where L(p) is a polynomial whose exact form will not concern 


tia tft () 
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us. When we put p=A in this identity, we obtain G(#)=AK(A), 
whence A=G(A)/K(4). Now K(A) can be found either by re- 
moving the factor p—A from H(p) and then putting =A, or by 
differentiating the identity H(p)=(/—A)K(p) and then putting 
p=A, which gives H’(2)=K(A). Hence we may state the follow- 
ing rule: 

To find the coefficient A of the term A/(p—A) when p—2 is a 

simple factor of the denominator, strike out the factor p—2 from 

the denominator and then put p=) in the expression which 
remains. Alternatively, differentiate the denominator with 
respect to p and then put p=A in the resulting expression 

G(p)/H'(p). 

Usually (but not always) the first version of this rule is more 
convenient to apply. 

When p—A is a multiple factor of H(p) (i.e. when m>1), we 
first put p—A=u and express G(p)/H() in terms of w; it will 
then take the form L(u)/w"M(u) where L(u) and M(u) are 
again polynomials and M(u) does not contain the factor # (i.e. 
M(o)+0). We now want the coefficients in the expansion 

Lu) _ Am | Ama Ay 
u™M(u) a” irr ea aa 

Multiplying up by w™M(u), we get 

L(t) =(Am+A matt... +-Aye"-2)M(u)-+u"N(u), (7) 
where N(u) is a polynomial whose exact form will not concern 
us. We can now find Ay, Am-y,.++, A, (in that order) by 
successively equating coefficients for u°, u4,..., «1, The 
procedure will perhaps be best understood by seeing it at work 
in a numerical example. 

Example 4: Find the terms corresponding to p—x in the partial 

i P°+3P 
Sractions 1” aya V(b +3) 
Putting p—-1 =u, we have 
P'+3p=(4+1)9+3(u+1) =4+6u+3u2 +0 ! 
and (P+1)(b+2) =(u+2)(u +3) =6-+5u-+u% / 
So we want 


446 +307 +109 (A, + Agu +Au)(6+5u+u*) +... 
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the dots representing terms in u* and «*. We can get the correct 
constant term on the right by choosing A, properly; then the 
correct term in « by choosing Az, and the correct term in u® by 
choosing A,. Terms beyond x? do not concern us, and in 

the term u* on the left has no effect on the values of A,, 4, and 
A. The arithmetical work can be performed mentally, but those 
readers who mistrust their skill in mental arithmetic may prefer 
to set out the work more systematically like a long division (but 
done ‘in reverse’, working in ascending powers of « and ignoring 
terms beyond u*). Thus 


20 


Re tee 
+ 
t 


Peet 


What this calculations shows is that 
4+6u toutsuten (5+ Sut Sut) io sutw) Pie 
The fractions which we want are now got by dividing, 
2 
Faria by u* and then replacing « by p—1; thus they are 


2/3 4/9 1/54 
@-1* Gy p= 


We may now state the following rule: 


To find the terms involving p— A in the partial fractions 
G(p)/H(p) when p—A is a Latebepont the Ssseetad tal 
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put p—2=u and write G(p)/H\(p) in the form L(u)/umM(u). 
Then find Ay-+Amu-+ ...-+Ayw"- from equation (7), 
divide by u™, and finally replace u by p—A again. Te Stee 
Ag+... +A,0"-}, the method used in Example 4 may be 
imitated, u™ and higher powers of u being ignored. 


Naturally one must be prepared to meet complex factors of the 
denominator H(p). In practice G(p) and H() will usually have 
real coefficients, so that with any complex factor p—A the con- 
jugate complex factor ~—A will also occur, and_with any 
fraction A/(p—A) the conjugate fraction A/(p— 4). In the 
evaluation of the solution we will therefore obtain terms like 


(454 +dzby) r= dots ActaaR (Act) 


and these can always be written as real combinations of 
cosines, sines and exponentials. It must be admitted that 
multiple complex factors may lead to rather unpleasant cal- 
culations, but fortunately they do not occur very often. 


2.4 Further examples. [We continue to use the notations x, 
%, ... for the values of x, Dx, . . . when t=0.] 


Example 5: (D*+4)4=0; 
initial conditions x9=3, %,=*,=%3=0. 


The operational form is 


(b*+4)e=3p1, 
so that 
eae 
? pra 
Now 


Pt +4=(6*—21)(p* +28) =(6 —1 —1)(-1 +4) (P +11) (+1 +7), 

so that the factors of the denominator are all simple. If p—A is 

any one of these factors, the coefficient of the corresponding term 

in the partial fractions for -3°°- may be found by differentiating 
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the denominator and then putting =A (see the rule on p. 49 


which gives 
sy] 3 
4P* Jpn, 4 
for all four choices of 4. Hence 
i: e O ALSRE Paaee <Aa 
is Pes toi pordi perai teen)! 


= SR (eH fi 


Bet cost +3 ert cos t=3 cosh £ cos #. 
Example 6: (D—1)*#=sin ¢; 
initial conditions %)=1, %,=2, %,=3. 
The operational form of sin ¢ can be found thus: 
i p ‘P(b +4) p 
=I(ett) = ee eas 
sin £=1(et) (pet )- (ot 1) seat 
Also (p—1)*=p'—3p*+-3p—1. So the operational form of the 
problem is 


(p- — Hie (at—set+30) #219*-30) +0) +55) 


PPP PED 
-(o- Pipher Senn Te 


whence 

*_Pi—p +p —ptr, 

Pp (p—1)(p2+1) * 
The complex factors p+i in the denominator give rise to two 
conjugate fractions, and the coefficient of 1/(p—i) is 
r I I+t 
G=i “@+2n21 4a 


The corresponding contribution to 7 is at plus the conjugate 
term, that is 


an{ "Fes t+isin o}= —2(cost— sin #). 
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For the terms corresponding to the multiple factor p—1, put 
p—1=u, This gives 
p'—Pr +P Prt _T+2u+4ur+y .. 
(p—1)(p?+3) wletautul) * 

dots denoting terms in «* and u‘ which we shall not need. Follow- 
ing the method used in Example 3, we obtain 

T+2u+4qut+ ... =¥(2+2u-+u%) 

“bdu(2-+-2u-put) 


+Sute +2u+u) +... 


Thus the required terms are obtained by dividing 5 + Sut Sua by 
u® and replacing u by p—1; the corresponding contribution to # is 


(; ms ease Be + fa say 
2 (p=1)* oot . 


Evaluating this, and adding on the pea found earlier, 
we get the solution 


x= —~(cos t—sin t) +(21242445)et, 
4 cee air 
Example 7: (D* +*)x =f (t) (w+ 0 being a real constant); arbitrary 
initial conditions x =%, Dx =2,. 
The operational form is 
(p* +0*)* = (*op* +p) +f (4). 


+e, I 
sat pat trail 


The required partial fractions, easily found, are 


a (“ 
rr Ge @) {n+4) x lay ; . ) 


Hence 


pttor p—io ptin 'p*+o* 2i0\p—io ptio 
so that we get 


snan{i(no—)th 4, f. {elolt—2) etal} /(a)ade 
0 
=#,008 ot-+ “sin at [isin eo(t—n)f(r)dr. 
o 
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Example 8: (D*+2D+5)a=e-%; 
initial conditions x)= —3, *,=2. 
The operational form is 


(p*+2p+5)= -(- sorta) +a10)+55,)u 


also p*-+-2p +5=(p +1—2i)(p+1+2i). Hence we have 
*__ pt + 10p +7 
Pp (P+2)(p 41-2) 41421)" 
In the partial fraction expansion 
3*-+10p +7 jae 8 B 
~@ FAG F1—2(p Hi Fai) p42 Pyia pyr ta 
we have 
4-[- por = 12-2047 
P*+2P+5 |pn-e 
ee 
(+2)@+1+21) | ie (+2) 5 


I 
4-445 tS" 


Hence 
pep Pg ee 
#-(4585+ exar cal pas =a) 
tien Cena | 
5 
fe" +5e-#R((—8-44)(cos 2t+isin 2)) 
spent cos 2f-+sin 2¢). 
Example 9: (D4+2D*+1)z=0; 
initial conditions %,=1, %,=0, ¥)=—1, %=0. 
‘The operational form is 


(t-+-ap*-+3)x = (2(pt+2p%) —1(p))1 = (p+ pr. 
We are lucky cca the factor (p?+1) can be cancelled, and 


1/2, ? 
tpt eet +p) “Hp5ten)* 


So the solution is 
x =}(c! +-e-#) =cos t, 
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Example 10: (D*+2D*+1)x=0; 
initial conditions %9=x,=%,=0, %=1. 
This time we have (p*+-2p"+1)*—pr, and eg The 
denominator is (p—i)* (p+i)*; the terms in p—i are found by 
putting p—i=—u, which gives 


i I I 
Pr wera wate) 
Now 1=(—}—fiu)(—4+4iu+ ...)+ ..., dots denoting terms 
in u*. Dividing —}—}iu by u* and replacing u by p—i, we find 
that 
5 rae Ee ues plus the conjugate terms. Hence 


(p—-#)? 
a= aR {—Helt—Jielty 
=—Hoost+4sint. 

We have deliberately set out the solutions to examples in 
such a way that only the formulae (1) and (2) from p. 42 need 
be used. More elaborate formulae could sometimes be used to 
shorten the work slightly, the most prs: being 


5d do 
Got eG G=« Eee g=esinot. (8) 
We do not advise the reader to make strenuous efforts to 
memorize (8), since it can be derived very simply by taking 
real and imaginary parts in 
P___panettat, 
p-—a—iw 


2.5 Simultaneous equations. These can also be treated by 
the operational method. First consider the simplest problem: 
to solve two first order equations 
(4D-+b,)2+(4D+4d)y=f), @) 
(a:D+-b,)x-+ (¢D-+4s)y=g(), 
subject to the initial conditions 
X=%q, Y=Yo, When =o. (zo) 
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The equivalent pair of equations involving Q is found by 
integrating equations (9), which gives 
(2%) +-5,0%-+-04(y—yo) +d, 09=O/(t), 
4(%—Zq) +b.0%-+Co(y—Y0)-+4,0y=Og(t). 
Now put Q=/-, and multiply by #, obtaining 
(ab+0))*+ (4p+4,)y=(*oty HG) PI +0, (x1) 
(Ax -+-b2)%+ (Cap +ds)y= (alta t+Iota)PI+8(6). 

This is the operational form of the problem; if the operational 
Keone of fo rite are keen, we of course insert these. 
e equations (11) are now by ordinary bra, i 
treated like a number, and the psgies ost pm ra mal 

way. 


Example 11: 
(D—1)x-2y =} 
—2*+-(D—1)y=t; 
initial conditions *)=2, yo=4.. 
The operational form is 
(b—1)y—2y =2p1-+2, 
—2x+(p—1)y=4p1 +h 
Solving for x, 
tem -ae=(0-n)(20-+5) +2(a045))x 


LS a at 
Pp PXp+1)(p—3) * 
After finding partial fractions in the usual way, we obtain 
1/3_3/9_ 1, 28/9) 
empl —2—-9 
#( ep pat p= ! 
ty eh ap a8, 
We could find y similarly, but it is slightly simpler to note that 
I 28 Sie 28 
2y=Dx—x—t= ——+e-t 4 — et p34 4g t_ “Spat f, 
ee iine S17) Soe or 
I 2 
whence y=——t—— + e- “at 
aie 
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In exceptional cases the above treatment may break down. 
Suppose that in equations (9) the terms containing derivatives, 
a,Dx+-¢,Dy and a,Dx+c,Dy, are proportional. It is then 
possible to obtain a single equation, involving x and y but not 
their derivatives. For example, if the equations are 

Dx-+-2x-+-Dy—3y: 
2Dx-+-x+2Dy+2y=0, 
then we can derive the single equation 
3x—8y=21. 

In such a case, the initial values %» and yp can certainly not be 
independently assigned (in the above example we must have 
3%)—8y=0), and even when x, and yp are suitably connected 
there may be no solution to the problem. We shall not treat 
such cases; they should not occur in physical applications, pro- 
vided that the physical problem has been correctly translated 
into mathematical terms. 

The treatment of m simultaneous first order equations follows 
lines which should now be obvious. Denote the unknowns by 
% +++» %p; let the equations read 


DY D+onde=h) (=, 2,-0.4m) (ra) 


onl 
and let the initial conditions be 
%,=a, when f=o (r=1,2,...,m). (13) 
Then the operational form of the problem reads 


a ® 
D Gnd +n) =>, OnePTHL ll) (P=1, ++ 4M), (14) 
wl ol 
and these equations are to be solved by ordinary algebra. The 
solution can then be evaluated in the usual way. It is assumed 
that the terms in (12) which involve derivatives are indepen- 
dent, i.e. that it is impossible to combine these equations so as 
to obtain a single equation free from derivatives. This assump- 
tion is equivalent? to the following one: 
1See P, M. Cohn, Linear Equations, Ch. V, § 12 (in this series). 
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The determinant formed from the coefficients ay, (r,8=1,...,) The coefficient of # on the left is 
is not zero, P*—5p* +5p*+5p —6 =(b +1)(b —1)(P—2)(P—3); 
and it ensures that the initial values a, . . . , % may be inde- on the right we have (0-119 -3)0*—4p-41), 
pendently assigned. ae ere re _P(b—1)(b—3)(b?—4P +1) 
More complicated systems of equations can also be dealt gah Rut abaya a P-4 i 
with. We content ourselves with one illustration: suppose we Hence we get 
have two second order equations in unknowns x and y. Let 7" 
%q) %) Yor J, Aenote the values of x, Dx, y, Dy, when t=o, and RAAT ola 2 -( B/S. 0/8 i 
let the first of the two equations read P (P+Np—2)P—4) (P41 pt p—4 
(aD? +-bD-+0)x-+(a'D*+b'D-+c'\y=f(d). pmters Loa 
Then the corresponding operational equation can be written 5 10 
down by an obvious extension of the rule (p. 45) for a single To find y, we have 
second order equation: q (p—1)y =p(p —4)a —(P*—4p +1)1 
on the left, replace D by p; on the right, add to f(t) the terms pt 
(valap*4+Dp) +a (ap) +-yola'p*+-0'P) +31(0'P)X- “0-494 Sip) 
The second equation is treated similarly, and the resulting pair p+ 
of operational equations is then solved by ordinary algebra. =(P*—4p +) GFn@—2) ? 
Example x2: (o*—4D\9=(D=1) Y_(P* 4p +P +2) | "Ga eh tt ide pa) 
—4D)x—(D—1)y=1, as a mS 
(odin PPP +1)(P=1)(P—=2) | \P PHI PtP 
initial conditions %)=%,=J9=1, J, =0. y=r—e-t+ 3et—2ett, 
The operational form of the problem is Once again, difficulties arise when the terms involving second 
(p*—4p)* —(p—1)y =(1(P?§ 4p) +1(p) -+1(—P) +1) derivatives are not independent. For two equations this means 
=(p*§—4p +1)1, simply that the second derivative terms are proportional; in 
p the example, these terms were D*x and Dy, and no trouble 
(P-+6)¥-+-(p*—P)y = (10 +1(p?—p) +0(P) i i occurred because these terms are independent. 
~(+58) 2.6 Justification of the method. We have not yet supplied a 
P-4, detailed proof of the facts that the operational method, applied 
Solving for x, either to an nt order equation (§ 2.2) or to a system of first order 
ODF Re a 6))x=(p2— “= equations (§ 2.5), gives the correct solution to initial value 
Uept—4p)(*—2) + (P—1\b +6) 4 = (PYF) ue problems. This gap will now be filled: the proofs will be given in 
+o—a(or4 5h terms of Q (more convenient than p in theoretical work). 
= Consider first the problem of solving the equations (12) on 


=(pe—4p*4+4p*— Moms 1) Pp. 57, subject to the initial conditions (13). If %,,..., %, are a 
re solution to this problem then by integrating (12) and using the fact 
58 E.D.E.—E* 59 
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that QDx,(¢) mae —#,(0) =2,(8) — rawr obtain 


Sentes ) + Doonan, 


coi 


Sent bine > tee +00) 
= rf 
=g,(t), say (r=1, . . .m), 
versely, if 4, . . ., % Satisfy (15) then (differentiating) they 
satisfy the differential equations (12). Also, putting ¢=o in (15), 
the terms involving Q vanish and we get 


(15) 


Danlel0) a) <0 (rt...) 
{ 


But we are assuming that the determinant of the a,, is non-zero, 
and therefore these equations imply that x,(0)—a,=0, i.e. that 
the initial conditions (13) are satisfied. Thus we have shown that 
the initial value problem (given by (12) and (13)) is equivalent to 
the set of equations (15). Note that when we put Q=p-1 in (15) 
and multiply through by p, we obtain precisely equations (14), i.e. 
the equations which are actually used in practical work (where 
p is more convenient than Q). 

‘We now show that the equations (15) have exactly one solution. 
Denote ap,-+-bysQ by ¢,,(Q), let A(Q) be the determinant whose (r, s)*® 
element is ¢,,(Q) and let C,,(Q) be the cofactor of ¢,,(Q). Then A(Q) 
is a polynomial of degree at most » in Q, say 

4(Q)=49+4,0+ . . +400", 

in which the constant term 4, is simply the determinant of the a,, 
and is therefore different from zero; also C,,(Q) is a polynomial in 
Q. (We are using determinants whose elements are operators; they 
can be defined and treated like numerical determinants because 
their elements all commute, so that multiplications can be per- 
formed in any order.) Now let %;, . . ., #, be any solution of (15). 
Multiply the rt” equation in (15) by C,(Q), and add; then the 
coefficient of x, on the left is 


. 
D-eiahentay ={412) wren sh 
r= 
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So only x, survives, and we must have 


A1aye1= D Col her Gat... m). 


‘The operator 4(Q) on the left has the form 4g(1 +4,Q+ -. -+4nQ”) 
with 4,+0; it therefore has an inverse (cf. § 1.5). It follows that 
we must have 


Srl) pe 
Aye) Gat (x6) 


A similar argument, which we omit, will show that (16) does give 
a solution to (15); this completes the proof that (15) has exactly 
one solution, given by (16). Of course all we have done in the above 
proof is to go through one of the methods of solving linear equa- 
tions with numerical coefficients and to check that the steps are 
still valid when the coefficients are operators of the form @ys-+bysQ 
the resulting formula (16) for the solution is ‘Cramer's rule’.! This 
justifies the second step in the operational method, which consists 
in solving (15) ‘by ordinary algebra’. The final step is to evaluate 
the solution (16), which involves putting rational functions of Q 
into partial fractions; this has already been justified in § 1.5. We 
have now established the following fact: 


the differential equations (12) have exactly one solution satisfying 
given initial conditions of the form (13), provided that the deter- 
minant formed from the coefficients ay, is different from zero, and 
this solution can be calculated by the operational method as set out 
‘in § 2.5 (in its practical form, using p rather than Q). 


We turn now to the initial value problem for a single mt» order 
equation, stated as (4) on p. 44. To justify the operational 
method of solving this problem, we shall replace it by an equiva- 
lent problem for a set of m first order equations, Let us take %, 
Dz, . . . D"-4x as new unknowns and denote them by q, %4, . +» 
My_y. Then we have Dity=t, Diy =U «+» Digna tina 04 th the 
differential equation for + (in (4)) becomes Dity_,+-@,4p—1+ + 
+Antto=f(t); also tt, . . ., tn, have the initial values %o, « 
Conversely if té, - ” 1 Sigse eatiafy. the’jprecoding equations aad 
initial conditions, then #=w, is a solution to problem (4). So an 

1 See P. M. Cohn, Linear Equations, Ch. V, §§ 12-13 (in this series). 
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initial value problem equivalent to (4) is given by 
Duip— ty =0, Dit, —tty =0, © « . Dity-g—ten—y =0, 
Dit p+ jtin—y + Aatln—at ++ +Antly=L(t); (17) 
Ho = Fo Hy =H, «+ +4 Mp1 =%y-y When t=0, (18) 
But this, when inessential differences in notation are disregarded, 
is precisely the type of problem stated in (12) and (13); the crucial 
determinant, formed from the coefficients of the derivative terms, 
is non-zero because its elements are 1 on the main diagonal and o 
elsewhere. So we may write down the operational form of the 
problem as in (15), namely 
tho Qty = oy thy — Qty Aap « «oy Una Qn =F nmay } (19) 
Una +41 Qten-1 + Og Qtin—a + +anQttg =An—1 +27 (t).. 
We may then solve these equations by ordinary algebra; only the 
solution for %(=*) is wanted, and this can’ be obtained by 
multiplying the equations (19) by 
T+QQ+ .. --+Gn1Q"-1, 
. Q+a,0?+ ... +anQ"-1, 
arrsegr, 
m1, 
and adding. Writing x for 1, this gives 
(1440 +4,0* +... +4,Q")x=Q"/(!) 
Holl +aQ+ ... +ay-1Q"-1) ete) 
+A(Q +H OQAF ». Gna") +... Qh 
and x can now be found because the operator 1-+-4,0-++ . . .-+@nQ” 
has an inverse. This justifies the operational method in its theoreti- 
cal version, using Q. It only remains to observe that when we 
replace Q by p-? in (20) and multiply through by ", we obtain 
precisely the operational form (5) of problem (4) which is used 
in the practical version (using p) of the operational method. We 
have therefore established that 
the initial value problem (4) has exactly one solution, and this can 
be calculated by the operational method as set out in § 2.2. 


2.7 The general solution of an n‘. order equation. The 
structure of the general solution of 
(D"+-4,D"-14 .., +a)x=f(b), (2x) 
becomes very clear when the operational method is applied. 
Whatever initial conditions are chosen, we obtain a solution 
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of the form 
C(p). I 
x=p——1+—— fi) (22) 
tat ag Bc 
where A(p)=/"-+-a,p""1+- ... 4+-@, and the polynomial C(p) 
depends on the initial conditions but does not depend on /(¢). 
Thus (22) displays the solution as the sum of two terms; the 


first term, p ae 1, satisfies the reduced equation 


("+ ... +4,)x=0 (23) 
with the same initial conditions as before, whilst the second 


term, ae {(), is the particular solution of (21) for which the 
initial values of x, Dx, . .. , D"~1w are all zero. Moreover we see 


that every solution of the reduced equation has the form 


ro when we evaluate this expression, after factorizing 
A(p) as 
(b—Ay™(b—Ag)™ « . . (P—Ax)™, 

we obtain a linear combination of 

ht, teh! emmteht: set, tmavtelit (24) 
Conversely, as we already showed in Ch. I (p. 21), any linear 
combination of the functions (24) is a solution of the reduced 
equation. We have therefore confirmed what was stated with- 
out proof on p. 22: 


the general solution of the reduced equation (23) is an arbitrary 
linear combination of the functions (24). 


Now suppose that we merely want to find the general solu- 
tion of (21). Then all we need, after writing down the comple- 
mentary function, is a particular solution, The operational 
method at once provides us with one particular solution, 


namely HH: We can do even better when the operational 


A(p) 
form F(#)r of f(¢) is known, so that we get the particular solution 
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THE OPERATIONAL METHOD 


ag because when we evaluate this (by partial fraction 
expansion) we may omit all terms which have already occurred 
in the complementary function. For example, a particular 
solution of (D*—x1)x=te' is 


Mee MBE Tyee : 
p§—1 (p—1)? “ (P*+-1)(b+1)(P—1)? 
but when we make the expansion 
pe A, A BG, 
(P*+1)(b+1)(P—1)? p-i pti’ ptr pr 
Cs Gs 
‘oO 


we need only the last two terms because the other terms will 
merely produce multiples of e*, e~*, e~* and e', which are 
already in the complementary function. The technique for 
partial fractions described in § 2.3 is very convenient here: it 
enables us to isolate a particular group of terms without having 
to find the other (unwanted) terms. The numerical work in- 
volved is perhaps no shorter than if one were to substitute 

trial solution, but has the advantage of following a familiar 
pattern when one has become used to the operational method. 


PROBLEMS FOR CHAPTER II 


[The values of x, Dx, D®x, . . . at t=0 will be denoted by x», 
*y %,...; D stands for d/dt.} 


1. Sone the partial fractions occurring in Examples 1-3 of 
2. 


§ 

2. (D?+6D-+9)x=1e-; xo=. 
3. (D?—D—2)x=e' sin 2t; x, 
4. (D*—5D*+5D*+-5D —6). 
5. (D*+-3D*—D—3)x=e~*; 
6. (D?+-2D-+-2)x=f(t) ; x=. 
7. (D*—2D+-1)a=e'; x, 

8. (D8+1)x=e-*; xy=1, 4,=—I, 2, 
9. (D?—5D+6)x=cos 3¢; %9=0, %=5. 
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3=2, X= %_=0. 
%_g=0. 


SOLUTIONS TO PROBLEMS 
10. (D?—3D?4-3D—1)x=54e"; %=I, 4=2, %=3- 
11. (D'4+-2D84-D*—2D—2) x=; x=, %y=—2, %_=0, X= 5. 
12. (D'—5D*+6D*+4D—8)x=0; xp=0, %=%yg=I, %y=3- 
13. (DY 49)x=/(0); m= m= M=H_—O- 
14. (D*-+-4D*+-4)x=1+0; By= m= Fy=%q=0. 
15. ease | Xp =I, Xy=O. 
16. +-2)x-+-Dy=0, yey 
Dit tab ayaa SmI i 
17. (D*—8)x-++-Dy=0, 
—6Dx+-(D*+2)y=0,, 
18. (D?+-5)x—y=4t, 2 
( ln ree Yomi—=0- 
19. D¥x+-2x—-2=D*y+2y—2=D%z—8x—y-+-82=0; 
X= =Yo' 1==0, 2y=1. Find z. 
20. (2D+3)x+- (oe 
( Fi NA Find x. 


Xp=JYo=HI, 4=—I, 1 =6. 


SOLUTIONS TO PROBLEMS 


[Arbitrary constants will be called A, B, C, sometimes with 
suffixes attached.) 


Chapter I (p. 25): 
1. Ce*—(x*+-2x+2). 2. xe. 
3. ¢cos x+} sin x+-Ce-™, 4. e~*(C+-log(1+-e™)). 
5. e~**(2 +-A cos x-+-B sin x). . 
6. fe-*(cos x— 2 sin x)-+-A cos x+-B sin x. 
7. 4 cos x-+Ae*+-Be*, 8. —}x*+-fx+-44+-Ac*+Be-™, 
9. (A+Bx+-4x)en*. . 
10. pyx(2 sin x—cos x)+-y'y(cos x—7 sin x). 4 
place R.H.S. by xe, substitute trial solution 
(Ax+B)e, then take the imaginary part.) 
M1, $e-*#—}e™. 12. }x(cos.x+-3 sin x). 


13. ie raerioe(B cos Viv4¢ sin %s), 


14. 2 cos 3x—$ sin 3x-+Ae~*+e**(B cos x-+-C sin x). 
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SOLUTIONS TO PROBLEMS 


par TI (p. 64): 


St, 


6. fie et) sin (t—z) f(r)dr. 7. e'(1—2t+-H). 


8. Hote 0(§ 0s (; via) fom cin(22)), 


9. at 8t_. 
abe abberstabe ah 008 a5 sn 39, 


41. 4(1—t) +-e~"(§ cos ¢-+4-sin t—2), 12. 4(e%#—e-4, 


13. if {cosh(¢—r)sin(¢—z) —cos(¢—z)sinh(t—z)}/(z)dr. 


14. 4(0°—1+-cos(t 
te irene ay 2 sin 2t), 


- t=aly(—Sel+ 176~t—2eWi"), y= aly (25et—17e-+Be-®), 


17. x=}e"+ fe-*+- (cos 2t—sin 21), 
=the" —e-*t+-3(cos 2t-+-sin 2). 
18, ae catty fom, 
=f? 2) G 
19. 2=y%5—} cos(t/2) +44 ea 


20. af" emf tetenae +] e~-(39(z)—2/(r))dr. 


2 e+ eet fot aoe an bores i par hve 
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Arbitary constants, 2, 17, 24 
Associative law, 31 


Commutative law, 31 
Complementary function, 5, 10 


Distributive law, 31 
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Initial conditions, 2, 17, 25 
Initial value problem, 27, 28 
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Operational form: 
of an initial value problem, 44, 56 
of f(), 45 
Operator: 
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polynomials in an, 32 
powers of an, 32 
product of two, 31 


Partial fractions, 37, 48 
Particular solution, 5, 10, 21 
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.>~G., EH. REUTER 


In this explanation of the simplest methods for, solving 
» differential equations with constant coefficients; Cl 

tajns a detailed description of standard mettiods for fin 
general solution and Chapter II gives an elemen 
for, finding the’ solution which satisfies spe 
ditions. 
Readers will find a careful description of practical 
of solution, and those who also wish to learn on what | 
matical principles these practical methods rest il fin 
"these Principles ‘are ‘discussed in an informal fe 
rigorous manner. « 
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